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I. INTRODUCTION 



Recently we initiated a program of calculating of all three-loop radiative-recoil correc- 
tions to hyperfine splitting (HFS) in muonium [1]. Three-loop radiative-recoil corrections 
are enhanced by the presence of the cube of the large logarithm of the electron-muon mass 
ratio [2]. All leading logarithm cubed and logarithm squared contributions of this order were 
calculated a long time ago [2-4] (sec also reviews in [1,5]). As the first step of our program 
we obtained in [1] previously unknown single-logarithmic and nonlogarithmic radiative-recoil 
corrections of order a^{Za){m/ M)Ep^ generated by graphs with two-loop polarization in- 
sertions (irreducible and reducible) in the two-photon exchange diagrams. As the next 
logical stage in implementing our program, wc present below the calcuulation of all single- 
logarithmic and nonlogarithmic three-loop radiative-recoil corrections generated by diagrams 
with one-loop electron and muon polarization insertions in the exchanged photons. There 
are four gauge invariant sets of such three-loop diagrams, and we calculate all their contri- 
butions. 

II. RADIATIVE-RECOIL CORRECTION OF ORDER a{Za){m/M)EF 

All four sets of diagrams considered below can be obtained from the two-photon ex- 
change diagrams with the radiative photons in the electron or muon lines by insertions of 
the one-loop electron or muon polarization operators. As was discussed in [1], it is suffi- 
cient to calculate contributions of these diagrams in the scattering approximation. In the 
calculations below we use the approach developed earlier for analytic calculation of the two- 
loop radiative-recoil corrections of orders a{Za){m/ M)Ep and {Z'^a){Za){m/ M)Ep in [6,7] 
(these corrections were also calculated numerically in [8]). To make this paper self-contained 
we first briefiy remind the reader of the main steps in calculation of the corrections induced 
by the radiative photon insertions in the electron line. 

The integral representation for the radiative corrections of order a{Za){m/ M)Ef gen- 
erated by the graphs with radiative insertions in the electron line in Fig. 1 has the form 
[9,10] 

SE'--" ^ aiZa)E, J^- [ f j -,J«^ (2) 



We define the Fermi energy as 

EF = ^Za-[-) chR^, (1) 

where m and M are the electron and muon masses, a is the fine structure constant, c is the velocity 
of light, h is the Planck constant, R^o is the Rydberg constant, and Z is the nucleus charge in 
terms of the electron charge (Z = 1 for muonium). The Fermi energy Ep does not include the 
muon anomalous magnetic moment which does not f actor ize in the case of recoil corrections, 
and should be considered on the same grounds as other corrections to hyperfine splitting. 
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FIG. 1. Electron-line radiative-recoil corrections 
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where the dimensionless integration momentum k is measured in units of the electron mass, 
the small parameter /j, is defined as half the ratio of the electron and muon masses {/j, = 
m/(2M)), the auxiliary functions of the Feynman parameters a{x, y) and h{x, y) are defined 
by the relationships 



X 



X 



(3) 



and explicit expressions for the coefficient functions q are collected in Table I. We used the 
Yennie gauge for the radiative photons in the derivation of the explicit expressions for these 
coefficient functions q. 

After the Wick rotation and transition to four-dimensional spherical coordinates, the 
expression in eq.(2) acquires the form 



5E 



•e-iine ^ a{Za)_m~ 

7r2 M ^ 87r/i 



1 rn roo 

— / d9 sin' 9 / dk' V(k,9) 

TU,'' Jo Jo 



(4) 



where 



k' 



k' + ii-' cos' 9 {k' + a')' + 4b'k'cos'9 ' 



V{k, 9) = dx £ dy{ (2 + cos^ 9) 



{c^ sm' 9 + c^k') (^)' + 2c^^] {k' + a') (5) 



-86c4(2 + cos^ 9) cos' 9^ + 126 cos^ 9(05^- cr) - 12ce cos' 9^(k' + a')}. 



We introduced derivatives with respect to a' in order to reduce the powers in the denomi- 
nators before integration over angles. 
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Next we separate the //-dependent and //-independent terms in the integrand with the 
help of the identity 



( + /(-2 cos2 e)[{k^ + a2)2 + 4 62 A;2 cos2 e 



(6) 



(A;2 + a2)2 - 4 /(2 62 ^4 



/(2A;2 + cos2^ (A;2 + a2)2 + Ah'^k'^cos^e 



1^' 



Ab'^k^ 



//2A;2 + cos2^ (A;2 + a2)2 + 462A;2cos2^ 



where, at the last step, we omitted the term proportional to /x^ in the denominator before 
the square bracket since, as can be shown, this term generates recoil corrections which are 
at least quadratic in the recoil factor /i. Then the integral for the radiative corrections in 
eq.(4) becomes a sum of /t-dependent and /{-independent integrals 
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(7) 
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H^k^ + cos^e {k^ + a^y + Ab^k^cos^e 



+ 5e: 



e—line 



which we will further call /i-integrals and c-integrals, respectively. 

It is convenient to write the angular integrals in eq.(7) in terms of the standard functions 
$1 defined as 



^n{k) 



TT/J,^ Jo 



d^sin^^cos^"^- 



,2\2 



(A;2 + //-2 cos2 e) [(A;2 + a2)2 + 462 A;2 cos2 9 



(8) 



TT Jo 



desin^ecos^'^e 



1 



462 A;2 



l^^k^ + cos^e (A;2 + a2)2 + 462A;2 cos2 ^ 



Explicit expressions for these functions and their properties are collected in Appendix A 
(see also [10]). 

The c-integrals with the functions ^n{k) do not contain any free parameters and generate 
only contributions linear in m/M to HFS. These integrals are pure numbers, which admit 
analytic calculation. The /i-integrals with the functions $f (/;;), ^n{k) parametrically depend 
on II and generate both nonrecoil and recoil contributions to HFS. Contributions of a fixed 
order in the small mass ratio (which are often enhanced by the large logarithms of the 
mass ratio) can be extracted from the /t-integrals with the help of an auxiliary parameter 
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a chosen such that it satisfies the ineqiiahty 1 <^ a <^ The parameter cr is used to 

separate the momentum integration into two regions, a region of small momenta < A; < a, 
and a region of large momenta a < k < oo. In the region of small momenta one uses the 
condition /ik <^ 1 to simplify the integrand, and in the region of large momenta the same 
goal is achieved with the help of the condition k^ 1. Note that for /c ~ cr both conditions 
on the integration momenta are valid simultaneously, so in the sum of the low-momenta 
and high-momenta integrals all cr-dependent terms cancel and one obtains a cr-independent 
result for the total momentum integral (for a more detailed discussion of this method, see, 
e.g., [11]). All contributions, up to and including the corrections quadratic in the small mass 
ratio m/M, were analytically calculated earlier in this framework [6,10] (see also [11,5]). 

III. DIAGRAMS WITH RADIATIVE PHOTONS IN THE ELECTRON LINE 
AND ELECTRON VACUUM POLARIZATION 

The general expression for the contribution to HFS arising from the diagrams in Fig. 2^ 
is obtained from the integral in eq. (4) by insertion in the integrand of the doubled one- loop 
electron polarization (a/n)k'^If.{k) 

2(^)em=2(^)k' tdv f'']-f'\ , (9) 

Vtt/ ^' Kit J Jo 4 + F(l-'y2)' 

where the additional multiplicity factor 2 corresponds to the fact that we can insert the 
vacuum polarization in each of the exchanged photons. 
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FIG. 2. Electron line and electron vacuum polarization 

Then all contributions to HFS generated by the diagrams in Fig. 2 are given by the 
integral 

SE^' ^ JL r de smH Hde v(k,e) (lo) 

7r3 M AiTfi^ Jo Jo ^ ^ ^ ^ 

1 v'{l-vV3) 

F + ^-2 cos2 9 (A;2 + a2)2 + cos^ 9 Jo 4 + k'^{l - v^) ' 

To obtain the radiative-recoil corrections of order a'^{Za){m/M)Ep contained in this integral 
we follow the route described in the previous section. First we write the integral in eq.(lO) 



^The graphs with the crossed exchanged photons are not shown explicitly in this figure and similar 
figures below. 
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as a sum of seven /i- and seven c-intcgrals. These integrals are collected in Table II and 
Table III, where the dimensionless contributions 5e to the energy shifts are defined as 5Ei — 
5ei a\Za){m/M)EF/7r^. 

The c-integrals in Table III automatically contain only contributions hnear in the recoil 
factor m/M, and can be immediately calculated numerically. The respective results are 
again presented in Table III, and the total contribution of all c-integrals is 

SET ^ K' ^'-^E. - 6.96182 (^)^^E.. dD 

The situation with the /^.-integrals is more complicated. Besides the recoil contributions 
linear in the small mass ratio m/M, they contain both nonrecoil contributions and the recoil 
contributions of higher order in the mass ratio. We would like to remove the already known 
nonrecoil contributions, to extract in analytic form all coefficients before the logarithmically 
enhanced terms linear in the recoil factor, and to calculate numerically the nonlogarthmic 
term which is linear in the recoil factor. To this end, we need to throw away the recoil 
contributions of higher orders which are also contained in the //-integrals and which are 
too small from the phenomenological point of view. If preserved, these higher order recoil 
contributions result only in the loss of accuracy in the numerical integrations. It is easy 
to see that the integrals with the function $o(^) generate only already known nonrecoil 
contributions and can be safely omitted for our present goals. We extract the contributions 
linear (and logarithmically enhanced) in the recoil factor from the //-integrals by separating 
the integration region with the help of an auxiliary parameter a as described in the previous 
section. All logarithms of the mass ratio originate from the high momentum parts of the 
/i-integrals, where we can use the high energy asymptotic expansion of the polarization 
operator. Thus we obtain in the analytic form all coefficients before the logarithms of the 
mass ratio. The radiative-recoil corrections of order a{Za){m/M)EF in eq.(lO), generated 
by the insertions of the radiative photons in the electron line, are linear in the large logarithm 
of the mass ratio. Hence the corrections of order a'^{Za){m/M)EF are quadratic in this 
large logarithm. The coefficient before the logarithm squared, which we obtain in this way, 
coincides with the one obtained earlier [3], and the analytically obtained coefficient before 
the single-logarithmic term, as well as the nonlogarithmic contribution, are new. The results 
of the calculation of the /(-integrals are collected in Tables IV and V. 

As an illustration of our methods let us obtain the first three entries in Table IV. We 
start with the integral in the first line in Table II. It contains a nonrecoil contribution 

5e-inon - recoU) = dx dy c, d^ j^^.m) ^ + (12) 



' dv- 



1 k 
^JtJo "^^X Jo "^^'W+^Jo ^^4 + F(l-^;2)- 

This integral, as well as all other nonrecoil contributions of order a'^{Za)EF, was calculated 
analytically in [12], and numerically in [13], and we will not consider it here. We write 
the recoil part of the integral in the first line in Table II as a sum of low-momentum and 
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high-momentum integrals^ 

]^ rl i-x f-oo 

2 Jo Jo Jo 



+ a2)3 



-I Jo 



t;^(l-t;V3) 
4 + A;2(l-i;2) 

(13) 



First we estimate the high momentum part 5e^^{' considering the leading term in the inte- 
grand in the high momentum hmit (more technical details on similar estimates can be found 
in [10]) 



1 r°° k'^ r 1 



z;^(l-t;V3) 
4 + k\l-v^) 
(14) 



°° dk^ 
■ 1^ 



2 u2 



In'k 



2$^(fc) - ^^,{k) - ^^{k) 



In^ a. 



We see that the high momentum contribution is suppressed like l/cr^. We have already 
extracted the recoil factor from 5e^^^ explicitly, so now we are looking only for such con- 
tributions which are not additionally suppressed. Hence, in the leading order in the recoil 
parameter we can omit the high- momentum contrbution to Se^i, and in our approximation 
the total contribution to the energy shift is given by the low-momentum integral 



t;^(l-t;V3) 
4 + k\l-v^)' 
(15) 



To extract the leading nonvanishing contribution to this integral at /x ^ we substitute 
in the integrand the leading terms in the small n expansion of the functions ^i{k) (see 
Appendix A) and obtain 



r ee< 



1 

4 + km 



(16) 



The momentum integral is ultravioletly convergent, so the magnitude of the integral with 
an infinite upper limit differs from the integral in eq.(16) only by inverse powers of a, which 
we omit anyway. Hence, the total /il contribution to the energy shift is given by the integral 



Here, as well as in similar cases below, we allow ourselves a slightly confusing notation, using the 
same symbol Se^^i both for the total contribution of the /x-integral and for its recoil part. We hope 
this will not lead to any misunderstanding since in this paper we are interested only in the recoil 
corrections. 
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(17) 



which is just the first entry in Table IV. 

Let us turn now to the next integrals in Table II. Each of the integrals (5e^| and (5e^| 
generates spurious logarithm cubed contributions which cancel in the sum of these terms. It 
is therefore convenient to rearrange these terms in such way that spurious logarithm cubed 
terms do not arise at all. To this end we write the coefficient functions C2 in the form 



C2 = 



4(1 - 2y) 

y{^-y) 



,t(-1 + x + 21nx) + 



4 



yf 



-{l-x){-l-2y/x) 



(18) 



+2(-4 + Ay/x) Inx + x{-l + x + 21nx)(4 - 5y + 2y^ 



= C2a + C2b- 



We have chosen C2a essentially as the singular part of the function C2 at y ~ and multiplied 
this singular part by the factor (1 — 2y)/{l — y) in order to simplify integration over y. 
Since (1 — 2y)dy / {y{l — y)) = —da^/a^ this makes integration over y trivial. The spurious 
logarithm cubed term originates only from the integral of the function C2a- On the other 
hand, the integrals over x and y with the function C26 remain finite even if we omit a? in 
the denominator (A;^ + a^)^, which leads to significant simplifications in the high-momentum 
part of the integral 5e^|. 

Following the same logic, we represent the coefficient functions C3 also as a sum of two 
functions and c^b 



y{^-y) {^-yy 



2(1 - y)(l - 6x - 2^2) + (1 - y)(8 + 2x) (19) 



-(l/a;)(26(l - x) + 2x{l - x) - 6(1 - x^)ylx - 6x(l - y) + 161nx) 



= C3a + C3&. 



Now the sum of the recoil contributions to HFS generated by the /x2 and jiZ integrals in 
Table II may be written in the form 



ee I r ee 
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Let us consider calculation of the (^e^2a+/i3a contribution 
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+ e 



ee< 



^ee> 



As mentioned above, calculation of the y integral becomes trivial after the change of variables 
y ^ a^. In the high- momentum integral we next expand over the inverse powers of A;^, and 
use the asymptotic expansions of the arising functions V;„^ (these functions were defined 
and their properties were discussed in [10]) 
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In the low-momentum part of the integral we again preserve only the leading terms in the 
small II expansion of the functions $f (A;) 



and integrate explicitly over y 
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(24) 



3 /'I 

— dx dk^' ' ' ' "'-^^ 

4 Jo io 



-1 + 6a: + 2a;2) 



lnA;2-ln(A;2 + a2) 



fc2 + af 



-4a;(-l + X + 21na;) 



-Ink^ + ln{k^ + al) + 



2aj 



k^ + aj 2 (A;2 + af )2 



1 , v\l - vVS) 
dv- 



4 + A;2(l - ^;2) j ' 

where = a;/(l — x). Note that we again extended the momentum integration to infinity 
in the finite integrals for the low-momentum contribution. Finally, summing the high- and 
low-momentum contributions we obtain the cr-independent result 
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which is convenient for numerical calculations. 

As our final example let us derive an expression for ^e^l^ in the third line in Table IV. 
We start with the general expression for this contribution from Table II 
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As we already mentioned above in this case we are free to omit in comparison with k'^ in 
the denominator in the high-momentum part of this integral 
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Deahng with the low-momentum part of the integral in eq.(26) we would like to extract 
analytically leading logarithms of a. This woTiId allow us to get rid of any trace of the 
parameter a in the final expression for this integral, making it much more suitable for further 
numerical calculations. The double integral over x and y, with the integrand including the 
function C2b, cannot be calculated analytically as easily as the integral with the function 
e^^2a+ij,3a above. In order to overcome this difficulty we use the identity 



k' 



(F + a2)3 F + 4 
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k' 



(k^ + a^y F + 4. 



(28) 



in the low-momentum part of the integral in eq.(26). Analytic calculation of the integral 
with the first term on the right hand side (RHS) is simple. On the other hand, all logarithms 
of a are supplied by this integral since the second term on the RHS in eq.(28) decreases at 
large k faster than the first term. Then we obtain the low-momentum contribution to the 
integral in eq.(26) in the form 
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Finally, the total expression for the fj,2b integral has the form 
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Calculation of the remaining /^-integrals goes along the similar lines, and does not require 
any additional new tricks. Calculating the nonlogarithmic contributions numerically and 
collecting all //-integrals from Table IV, we obtain the total contribution of all //-integrals 



5, 2^ 22, M 

-In^ h — In h4. 45606 1 . 
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(31) 



The sum of /t- and c-integrals in eq.(31) and eq.(ll) gives all radiative-recoil corrections of 
order a^{Za){m/M)EF generated by the diagrams in Fig. 2 



^ln2^ + ^ln^ + 11.41788 (3). 
2 m 3 m 



(32) 



IV. DIAGRAMS WITH RADIATIVE PHOTONS IN THE ELECTRON LINE 
AND MUON VACUUM POLARIZATION 



Let us consider now the diagrams in Fig. 3. The only difference between these diagrams 
and the diagrams in Fig. 2 from the previous section is that they contain the muon vacuum 
polarization insertion^ 



2 g) .%(.) ^ 2 g) 



1 ^, .^(1-.V3) 



(33) 



instead of the electron vacuum polarization, and the respective contribution to HFS has the 
form (compare eq.(lO)) 



TT^ M 47r/t 
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^We ascribe an extra factor Z to each photon emission by the heavy line, but being somewhat 
inconsequential, do not write any Z factor in the muon vacuum polarization. 
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FIG. 3. Electron line and muon vacuum polarization 



dv- 



[i-vys) 



+ 11-^ cos2 e (A;2 + a2)2 + 462A;2 cos2 9 Jo /i-^ + k^{l - v^) ' 

We transform the integrand in the sum of fi- and c-integrals exactly as in the previous section, 
the only difference being that each integrand now contains as a factor the muon vacuum po- 
larization insertion from eq. (33) instead of the electron vacuum polarization insertion. This 
immediately leads to significant simplification of further calculations in comparison to the 
case of the electron polarization. The muon polarization insertion in eq.(34) is suppressed 
as A;2|i2 integration momenta much smaller than the muon mass. But the characteris- 
tic integration momenta in the c-integrals are of order one in our dimensionless units, and 
hence now all c-integralsdo not generate corrections linear in the recoil factor m/M. Low- 
momenta parts of the /x- integrals, where integration goes over momenta k < a 1/ /i, 
are also suppressed as a^n^. Only the high-momenta parts of the //-integrals, where inte- 
gration effectively goes over the momenta comparable to generate contributions linear 
in m/M = 2/i. Such contributions are present only in the //-integrals connected with the 
coefficient functions C2, C3, cg, and C7. Note that the integrals with the coefficient functions 
C2 and C3 contain the function $o(^) — ^/il^^) integrand. In the case of the electron 

vacuum polarization (and absence of any polarization at all) the integrals with this function 
generated nonrecoil corrections only. This essentially happened because the characteristic 
integration momenta were of order one, and in this region the contribution of the function 
l/{fik) is obviously enhanced as In the present case, characteristic integration momenta 
are about 1/ (/i/c), and the integrals with the function $o(^") generate recoil corrections on a 
par with the integrals with the other functions $(A;). 

Integration over the Feynman parameters ,t; and y in the high-momenta integrals simpli- 
fies and may be performed analytically, and we obtain all one-dimensional momenta integrals 
which generate radiative-recoil corrections in the case of muon polarization insertions 
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1 , v\l-vV2.) 
dv 



/i 2 + k'^[\ — v"^) ' 



r 



dk^ 



4, , 5 

-Ink 1 

3 3. 



27r2 89 \ 
■^+18) 



2%{k) + 2$(;(A;) + ^tik) 



(36) 



1 v\l-v^/2,) 

''/(-2 + F(l-v2)' 



12 



TT 



21 

T 



)/; 



1 v\l - 7773) 



(37) 



1 ^. ^;7l-t^V3) 



/x-2 + A;2(l-^;2)- 



(38) 



The sum of all these contributions has the form 
or explicitly 



+ 3m)\l"dv (39) 



jj, -\- k'^{l — ^2) 



/•oo f 7 



(40) 



+2>[-llk^Jl + lJ?k^ + /^2A;2 + I ^ 



1 t;2(l-t;2/3) 
//-2 + A;2(l-^;2)' 



We can extend the momentum integration region in this integral to zero, since contributions 
from small momenta are additionally suppressed by powers of cr/i. It is also natural to rescale 
the integration momentum k iik. Then the expression for the radiative-recoil corrections 
has the form 



5eT 



no 



dka 



(41) 



+Qk(-k\flT¥ + + 2 



1 ^2(1 - ^;2/3) 

l + A;2(l-^;2)- 



The integral in the first line of this equation is proportional to the integral for the 
muon polarization contribution of order a{Za){m/M)Ep in [10], and the contribution of 
the integral in the second line may be calculated in the same way. Finally we obtain 



57r2 1_ 
12 ^ 18 



- 4. 0567796 . . . 



(42) 



V. DIAGRAMS WITH RADIATIVE PHOTONS IN THE MUON LINE AND 

MUON VACUUM POLARIZATION 

The expression in eq.(2), for the radiative corrections generated by the diagrams with 
the radiative photon insertions in the electron line Fig. 1, was obtained without expansion 
in the mass ratio. Hence, after the substitutions m M and a — > Z'^cn it goes into the 
expression for the corrections generated by the diagrams in Fig. 4^ 



^We use the convention that the heavy line charge is Ze. 
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+ 2 



+ 



diagrams with 
+ I crossed exchanged 
photon hnes 



FIG. 4. Muon-line radiative-recoil corrections 



-line 



[Z^a){Za) m 



, , Ep \ ( dx [ dy 
M AJo Jo 



d'^k 



(43) 



+ 



A;2 + 4//A;o + iO k"^ - A/iko + iO 



2k2) 



cik^ + C2(A;2)2 



+ 



c^k^ + C42ko 



( - A;2 + 26A;o + 0^)3 ( - A;^ + 2&A;o + )2 



C5/c^ + C6k'^2ko 
( - A;2 + 26A;o + )2 



+ 



A;2 + 2&A;o + 



The dimensionless integration momentum here is measured in muon mass units, and appar- 
ent incomplete symmetry with the expression in eq.(2) is due to asymmetry of the parameter 
H = m/(2M). 

Unhke the case of the radiative photon insertions in the electron line in Section 11, the 
integral in eq.(43) does not generate nonrecoil contributions [7,14]. This happens because the 
radiative photon insertions in the muon line suppress the contribution from the integration 
momenta of order of the electron mass which were responsible for the nonrecoil correction 
by an extra factor /i^. Hence, for calculation of the leading order recoil corrections we can 
safely allow // — > in eq.(43) 



SE^- 



TT^ M^^ 2 70 ^"^io in'(k' + iOr 



(44) 



3/uq 2k 



2\2 



k^ + 2bko + a2 )3 



+ 



k^ + 2bko + a? Y 



3ko 



csA:^ + CQk'^2ko 



+ 



k^ + 2bko + a? 



+ 2hko + a? )2 

After the Wick rotation and transition to the four-dimensional spherical coordinates we have 

1 



-line 



{Z'^a){Za) m ^ 
7r2 M ^27r 



1 fK fOQ 

- / d^sin^^ / dk^V{k,9) 
In Jo Jo 



(F + a2)2 + 462A;2cos2^' 



(45) 
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where the differential operator T>{k,6) was defined in eq.(5). There are no /i-dependent 
terms in eq.(45), and this expression for the energy shift is similar to the c-intcgrals in 
eq.(7). It can be formally obtained from the expression in eq.(7) by the substitution 



+ a2)2 



All corrections to HFS described by the integral in eq.(45) were analytically calculated 
in [7]. Contributions to the energy shift generated by the diagrams in Fig. 5 are obtained 
from the expression in eq.(45) by insertion in the integrand of the doubled muon vacuum 
polarization 

.i^)em..i^)ej:,.^^^. ,4a) 




FIG. 5. Muon line and muon vacuum polarization 

This expression is apparently different from the respective expression in eq.(33) because now 
our dimensionless momenta are measured in terms of the muon mass. Finally, the integral 
for the contribution to HFS generated by the diagrams in Fig. 5 has the form 

= -{Z--) {Za) rn I T dO s.n^ 6 H dk^ V(k,e) (47) 
TT'^ M t: Jo Jo 



(A;2 + a2)2 + 462A;2cos2^ Jo A + k^{l-v^) VV ' 7 M ^' 

All seven contributions are collected in Table VI, and the total contribution to HFS generated 
by the diagrams in Fig. 5 is 

56^" = -1. 80176 (2). (48) 



VI. DIAGRAMS WITH RADIATIVE PHOTONS IN THE MUON LINE AND 
ELECTRON VACUUM POLARIZATION 

Let us turn now to the diagrams in Fig. 6. The only difference between these diagrams 
and the diagrams in Fig. 5 is that now we have electron polarization insertions instead of 
the muon polarization insertions. The respective analytic expression is obtained from the 
one in eq. (47) by substitution of the electron polarization instead of the muon polarization 



15 



2k 



FIG. 6. Muon line and electron vacuum polarization 

.,2n _ „,2/Q\ _ ^1 



io 4 + - v^) ^ io 4(m/M)2 + - v'') ' ^ ^ 



The typical integration momenta in eq.(45) are of order of the muon mass and, since we are 
calculating to linear accuracy in m/M, it is sufficient to substitute instead of the electron 
polarization operator the leading terms in its expansion over m/M 

2A;2 / dv — — -V^ r- a; - In A;^ + - In . (50 

Jo 4(m/M)2 + A;2(l --t;2) 3 3 m 9 ^ ' 

Then the total contribution to HFS generated by the diagrams in Fig. 6 is given by the 
integral 

5E^^ = ^fd^sin^^ rdeVik,e) (51) 

77^ M 277 Jo Jo 

1 /2, ,2 4, M 10 

- In /c^ + - In — - — 



(A;2 + a2)2 + 462A;2cos2^ \3 3 m 9 

where the contribution ^e^*^ = J2i ^ct corresponds to the first term in the last brackets, and 
the contribution 5e^^'^ corresponds to the second and third terms in the same brackets. We 
calculate 5e(f^ numerically and collect all seven contributions to this integral in Table VII. 
The sum of all these contributions is 

5ef = 12. 94447 (4). (52) 

The integral for 5ef is proportional to the integral for the radiative-recoil corrections 
generated by the radiative photon insertions in the muon line, which is known analytically 
[7,14]. Hence, we can immediately put down an analytic result for 5e^^ 

M 10\ / 9 „ 2. ^ . 39 



Then the total contribution to HFS generated by the diagrams in Fig. 6 is equal to 

13 \ , M 



(^6C(3) - 47rMn2 + y ) In— + 24.32115 (4). (54) 
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VII. DISCUSSION OF RESULTS 



Collecting the results in eq.(32), eq.(42), eq.(48), and eq.(54) we obtain the three- loop 
single-logarithmic and nonlogarithmic corrections generated by the one-loop electron and 
muon polarization insertions in the exchanged photons 



where 



SE = Se' 



5e' 



M 



Ep + 



M 



E, 



5, oM 22, M 

- In' — -f- — In — + 7.36110 3 , 

1 m 6 m 



(55) 



(56) 



5e^^( 6C(3) - 47r'ln2 + — ) In— + 22.51939 (5). 

\ 2 J m 



(57) 



If we recall that in the real muonium Z = 1 and that the logarithm squared term was 
already calculated earlier [3] then the new three-loop single-logarithmic and nonlogarithmic 
corrections obtained above may be written as 



5E 



f 



6 C(3) - 47r'ln2 + 



83 



In — + 29.88049 (6) 

m 



m 



TT 



3M 



E^ 



F- 



(58) 



Combining this result with the result of our earlier paper [1] we obtain all three-loop single- 
logarithmic and nonlogarithmic corrections generated by the electron and muon polarization 
insertions in the exchanged photons 



5Etot 

or, numerically 



47r'ln2 + 



67 
12 



In — + 9.59318 (6) 

m 



m 



TT 



3M 



E, 



SE = -0. 028 8kHz. 



(59) 



(60) 



This correction has the same scale as some other corrections to hyperfine splitting in 
muonium calculated recently [10,15,16]. All these results improve the accuracy of the theory 
of hyperfine splitting and affect the value of the electron-muon mass ratio derived from the 
experimental data [17] on hyperfine splitting (see, e.g., reviews in [5,18]). We postpone 
discussion of the phenomenological implications of the result above until the completion of 
the calculations of the remaining three-loop radiative recoil corrections. 
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APPENDIX A: STANDARD AUXILIARY FUNCTIONS 



Calculations described in this paper are greatly facilitated by the use of the auxiliary 
functions $n(^) (n = 0, 1, 2, 3) defined by the relationship 



1 r 



TT/J,^ Jo 



d9 sin^^cos^"^ 



( A;2 + cos2 e)[{k^ + a2)2 + 4 62 ^^2 cos2 e 



(Al) 



The integral over angles may be explicitly calculated, and the result of the integration is 
conveniently written as a sum 



where 



(A2) 



(A3) 



(A4) 



^ 16' 



2' 



- Cc{CcW{Cc) - I) - i 



(A5) 

(A6) 

(A7) 
(A8) 
(A9) 

(AlO) 



The standard function W{$,) has the form 



1 1 



1 

16' 



(All) 
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W{i) = yi+l - 1- (A12) 

and 

One may easily obtain asymptotic expressions for the function W{^) 

^mWiO^^, (A14) 

lim W{0 - ^. 

High- and low- momentum asymptotic expressions for the functions may also be easily 

calculated. Let us cite low-momentum expansions, which were used in the main text for 
calculation of the contributions to the hyperfine splitting of relative order /x^ 



^i(^) ^ - - l^k + (nkf, (A16) 
2 



« 1 - M!. (A17) 
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TABLE I. Coefficients in the Electron-Line Factor 



Cl 


16 


(1 — x)(a; — 3?/) — 2?/ In a; 




C2 


4 


— {l — x){x — y — 2y'^/x) + 2(x - % + 4y^/x) Inx 


2/(1-2/)'' 


C3 


1 


1 - 6a; - 2a;2 - (?//x)(26 - Qy/x - 37a; - 2a;2 + 12xy + 161nx) 


2/(1-2/)^ 


C4 




2x — Ax^ - 5y + 7xy ^ 


y(i-2/)^ 


C5 




' 6a; - 3a;2 - 8y + 2xy ) 


2/(1-2/)^ 


C6 






C7 
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TABLE II. |Lt-integrals for the diagrams in Fig. 2 



"Si 



1 Jo dx Jo dyci dk"^ jwte^ 



^2 \ Jo dx Jo' dyc2 /o~ dPj^^ [2n{k) + 2$(;(fe) + ^1{k)] J^' d^gg^ 



■ k Jo dx Jo dyc3 J^ dk^ JWT^ 



2$S(fc) + 2${;(fc) + m)] Jo dv t^^2) 



Aj^dxJSdyhc,J^dk^j^J^[2^'i{k) + ^^^{k)\ J^dv^^^^ 



"1 ^„. ^^(l-^V3) 



6 J^ dx J,^ dybc, dk^j^^^^^ik) J,' dv^^^^^ 



3 J^ dx Jo- dye, Jo-' de^^^>l{k) Jo' d^gfe^ 



-3 J^ dx Jo^ dybcr /o~ dk^j^^^,{k) Jo' 
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TABLE III. c-integrals for the diagrams in Fig. 2 



Kl -\ lo dx !S dye, /o~ dk^ (ai.) 'f^^ de sin^ 9 

(2 - cos^ 9 - cos^ 9) (,.^„.).f 4^,. , /o^ dv^^0-^ -0.916881 (3) 



Kl -\ lo dx /o^ dyc2 /o°^ (ai^) 'pf^^ lo dO sin^ 6 

(2 + cos^ g) (,.^„./f4k..os^, /o^ ^-iS^ 1-76474 (1) 



-i dx !^ dye, dk^^T^l £ d9 sin' 9 

(2 + cos^ g) (,.+„.)g4t'.^cos^. /o rf^lrfep^ 5.93600 (5) 



(2 + cos^ g) (fe.+„.).f4gfe^cos^g /o dv :l%~(f_% 0.0312952 (2) 



<56l - 3 dx dybc, /o~ dk' ^ 1 d9 sin^ ^ 

^ (FwSra /o ^-3^ 0.0003483 (2) 



Setl 3 dx /(f dye, j'^ dk' ^ d9 sin' 9 

-OS' ^ oEwIIp^ Io' dv^^^ 0.0543901 (2) 



3 dx /(f 6cv /o°^ dfc^ i jjr ^.^2 ^ 

^o^' ^ (FWOT^ /o ^-S?^ 0.0919342 (3) 
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TABLE IV. Simplified //-integrals for the diagrams in Fig. 2 



16 Jo Jo "y<^l Jo "'^ Jo "-^ 4+k'^{l-v'^) 



-2.105912 (2) 



3 3 "T 1 Q Q "T IS 1« T" 



m ^ 18 18 



^ - 1 ) In^ ^ 



+! Jo Io° dk"^ <{ (-1 + 6a; + 2^^ 



lnA;2-ln (A;^ + af) 



- 4x(-l -l-ar-F 21na;) 



4_ZEI _ 
18 18 



2a- 



_ 1 a\ 



Jo '^^ A+k\l-v^) 



+0.4698395 (7) 



3 ^ 18 i m ^ 



9 ^ 27 j 1" m 



^ In^ 
3 ^ 18 j 1" m 



- ^ (307r2 - 271) (8 In 2 - 3) - f dx dyc^i, 



Jo dVT^ ^ 



4+fe2(l-!;2) 



27 ; m 



_57rf , 27l7r2 
18 108 



-0.797177 (4) 



47if _ 89 \ 1^2 M I / _ 356 \ M , 2n^ _ 89-11^ 

3 9 i TO \ 9 27 j 1' m 9 54 



4^ _ f ] ln2 M 
3 9 / TO 



+ 2^ (127r2 - 89) (8 In 2 - 3) + I dx dyc^b 



16t7^ _ 356 \ i„ M 
9 27 ) TO 



lo^dk' 



Jo "'^4+fc2(l-i)2) 



27r* 897r^ 



54 



5.398568 (7) 





2 /o '^a; /o dyhci /o dA;2 /q dv ^_^^.,^fJJ^ 


-0.3738824 (5) 


^^/i5 


3 /o /o dybc^ /o dk'^y^j^ /q ^t^4+i^(i_(^) 


0.127357 (3) 
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TABLE V. Simplified /Lt-integrals for the diagrams in Fig. 2 (continuation) 



3 2 ' ' Q "r Q I m ^ "T , 



m ' 18 36 3 



^ - I In^ M 
3 2 1 m 



+ (27r2 - 21) (8 In 2 - 3) + I /o dyce 



3 rl J rx 



fe4 _ fc2 



rl , i;2(i_^2y3-) 

Jo ^'"4+F(T^ 



I TT OTT 14 

"•"18 36 3 



-0.1944535 (6) 



3 2/ m V9 3/ m 18 36 3 



+ I ) In^ ^ 



- ^ (27r2 - 15) (8 In 2 - 3) - I /d /(f dy6c7 



fc4 



fc2 

F+4 



Jo «^'4+feii(l-j;^) 



III _ £i 4_ 10 
18 36 3 



-0.8472862 (4) 
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TABLE VI. c-integrals for the diagrams in Fig. 5 



<5er /o' dx dye, J,°° ( ^) ' i J,^ dO sin^ 9 

(2 - cos^ 9 - cos' e) ^^^^^^^^, dv^^i^ 0.297870 (1) 



<56f dx dyc2 /o~ de ( ^) ' i J,^ dO sin2 9 

(2 + cos^ ^) /o rf-iS^ -0-759522 (5) 



5e^/ 2 dx dye, dk^^^ £ d9 sin^ 9 

(2 + cos^ 0) ^,.^Xu:^^^e Io -1-11839 (2) 



5ef -8 /(I dy bc4 dP-^^^ d9 sin^ 9 cos^ 9 

(2 + cos^ g) (,.+..).^%fe.eos^. /o -0.0484183 (2) 



<56f 12 dx dybc, /o°° de ^ 1 /J^ d0 sin^ 9 

^os' ^Mr4lb'^k'^.o^e Jo dv^^l^^ -0.0027928 (4) 



.56^ - 12 dx J,^ dye /~ ^ i d0 sin2 9 

-s^ ^ /o ^-g^ -0-0251636 (1) 



5ef - 12 dx /(f dy6c7 /o°° dfe^ i £ d9 sin^ ^ 

(^WTW^^^e lo rf^lS^ -0.1453343 (4) 
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TABLE VII. c-integrals for the diagrams in Fig. 6 



\ J,' dx !^ dye, !^ dkH^]^- de sin^ e 



(2 



cos2 e - cos^ 9) In fe2 



15. 49349 (3) 



iJldxJ^dycJ^dkHD ^I,^d9sm'' 



(2 + co^e) ^.^Jp^^l,^^, h.k 



-1.235177 (4) 



I dx £ dyes lo^ dk' ^ i £ de sin2 9 



oo ,,2 9 1 . 



(2 + cos^ 9) i^T+^^r^^wW^^ 1^ ^"^ 



26. 74813 (3) 



- 1 !^ dx !^ dy he, /o~ de ^ i !^ d9 sin^ 9 cos' 

(2 + COS^ 9) ^kHaW^b^k-^oo.H 1^ 



1.304129 (5) 



4 dx dy6c5 /o°° dA;^ ^ i d0 sin^ 9 



0.411544 (4) 



-4 dx dyc6 ^ i /J^ sin2 9 



-0.0503852 (6) 



-4 /(J dx dybej dfe^ i de sin^ 



1.259721 (4) 
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